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Measurements of spatial diffusion of cold atoms in optical lattices iiave 
Q revealed anomalous super-diffusion, which is controlled by the depth of 

I , the optical lattice. We use first passage time statistics to derive the 

' diffusion front of the atoms. In particular, the distributions of areas 

swept under the first passage curve till its first arrival, and of areas under 
the Bessel excursion are shown to be powerful tools in the analysis of 
the atomic cloud. A rather general relation between first passage time 
statistics and diffusivity is discussed, showing that first passage time 
analysis is a useful tool in the calculation of transport coefficients. A 
^ brief introduction to the semi-classical description of Sisyphus cooling is 

iJ provided which yields a rich phase diagram for the dynamics. 

o 

1. Introduction 

in 

Pick's second law predicts how diffusion causes concentration to change 
in time. Instead of concentration we will use a probability density in one 
dimension P{x, t) and then Pick's law reads 



o 

o 



^ dP _ d^p 

where K2 is the diffusion coefScient. For an open system, and particles 
initially localized at the origin, the solution is a spreading Gaussian packet, 
with a width proportional to time: 

(a;2) = 2K2t. (2) 

Such normal diffusion is ubiquitous, which is hardly surprising since the 
diffusion is merely a sum of random displacements with zero mean. Hence, 
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following the central limit theorem we expect to find in experiment a Gaus- 
sian distribution for the particle position. The Gaussian universality of 
diffusion processes leaves us with one difficult task, and that is to compute 
the diffusion constant if 2 ■ This constant is non- universal, unlike the nearly 
universal shape of the Gaussian probability packet. Still, non-equilibrium 
statistical mechanics gives us a recipe (in principle) for the calculation of 
K2 which is the famous Green-Kubo formula 



where m is the mass of the particle. This relation between the momentum 
correlation function and the diffusivity is extensively used in the context of 
hard and soft condensed matter. The correlation function and the underly- 
ing process itself is assumed to be stationary so {p{t+T)p{t)) is independent 
of t. The Green-Kubo formula is important since once we obtain the dif- 
fusivity of a system we also know its mobility. This is made possible via 
linear response theory and the well-known Einstein relation which connects 
the diffusivity with the mobility via the thermal scale ki,T where T is tem- 
perature. 

All this was established a long time ago, so it is not surprising that physi- 
cists have turned their focus to problems that violate this normal behavior. 
One example is the diffusion and transport of atoms in optical molasses 
driven by counter propagating laser beams.l^ This system has attracted 
considerable attention, both for applied purposes like cooling and control 
of atoms and from a more fundamental point of view, since it exhibits un- 
usual friction and in some cases large deviations from ordinary statistical 
mechanics. Indeed laser cooling is today the tool of choice in many labora- 
tories which investigate low temperature physics. So basic understanding 
of statistical properties of laser cooled atoms is a timely subject. 

Within the semi-classical picture of Sisyphus cooling, the atoms are 
subject to a velocity dependent friction force, F(v), defined more precisely 
below. At low velocities, F{v) cx — u, which is of course normal in the 
sense that it mimics Stoke's friction for a macroscopic Brownian particle in 
a viscous medium at room temperatures. Such a friction force is intuitive 
since it indicates that the faster the particle moves the more energy it 
dissipates to the surrounding bath. The large v behavior of the atomic 
friction on the other hand is counterintuitive: it decreases with velocity, 
F(v) ^ — thus when t; ^ 00 the system becomes frictionless. Such a 
system can be called asymptotically dissipationless, and it exhibits physical 
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behaviors very different than standard systems such as non-dissipative ones 
and systems with kinetic dissipation which does not vanish at high speeds. 

This behavior of the friction force is due to the fact that the laser cooling 
is not effective beyond some critical velocity. Briefly, the effect is caused 
since a fast particle cannot distinguish the uphills from the downhills in 
the spatially periodic optical lattice; hence the Sisyphus mechanism which 
favors the loss of energy when the particle is close to the maximum of the 
periodic lattice is not effective (see Fig. [l] and caption) . Since the friction 
is vanishingly small for large velocities the equilibrium distribution of the 
velocities can be shown to obey power law statistics (see details below) , and 
so various average quantities are dominated by fast particles (compared for 
example with the Maxwellian distribution which is Gaussian). More gen- 
erally, the kinetics and equilibrium properties of atoms under laser cooling 
are far from standard: anomalous diffusion, non Gibbsian states, etc. can 
be found not only in Sysiphus cooling but also in other cooling approaches 
like sub-recoil laser coolingji^J and in Doppler cooling as well. 

Indeed in the closing paragraph of Cohen- Tannoudji's Nobel lectur^ 
he writes "It is clear finally that all the developments which have occurred 
in the field of laser cooling and trapping are strengthening the connections 
which can be established between atomic physics and other branches of 
physics such as condensed matter or statistical physics. The use of Levy 
statistics for analyzing sub-recoil cooling is an example of such a fruitful 
dialogue" . Here such a dialogue is extended to the popular Sysiphus cooling 
mechanism. Levy statistics roughly refers to power law distributions, which 
violate the conditions leading to the Gaussian central limit theorem. An 
example of this is found in a recent experiment performed at the Weizmann 
Institute. 

1.1. Spatial diffusion of cold atoms 

There, Sagi, et alW measured the diffusionSisyphus of ultra-cold ^^Rb 
atoms in a one dimensional optical lattice. Starting with a very narrow 
atomic cloud they recorded the time evolution of the density of the par- 
ticles, here denoted P(x,t) (normalized to unity). Their work employed 
the well-known Sisyphus cooling scheme.^^ As predicted theoretically by 
Marksteiner, et al.,'^the diffusion of the atoms was not Gaussian, so that 
the assumption that the diffusion process obeys the standard central limit 
theorem is not valid in this case. We recently determined the precise nature 
of the non-equilibrium spreading of the atoms, in particular the dynamical 
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phase diagram of the various different types of behaviors exhibited as the 
depth of the optical potential is varied, at least within the semiclassical 
approach. 
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Fig. 1. Sisyphus was punished by Greek gods to push a heavy stone uphill just to have 
it roll back to its starting point. In Sisyphus cooling, an atom will go up a potential hill 
created by laser fields, and then falls down to the local ground state, thus loosing energy. 
Such a scheme is made possible only when transitions from the potential maximum to 
minimum are statistically preferred, which in turn is made possible by clever quantum 
mechanics, and polarization fields (which control the transitions) which are periodic 
function of distance specially designed to correspond to the spatial modulation of the 
energy surfaces (the period is determined by the wavelength of the laser field, A^). 
However, if Sisyphus rans like a bolt, over the hills and the valleys he can fool the gods. 
Namely Greek gods cannot push him down from the top of the hills, since they need 
some finite time to implement his punishment. This corresponds to a fast atom whose 
transitions are not synchronized with it being on maximum of energy surface. So for fast 
particles, the Sisyphus friction becomes small and decreases, according to the detailed 
calculations, like F{v) oc —1/v. For more on Sisyphus cooling see.d 



In Ref. f4l the anomalous diffusion data was compared to the solutions 
of the fractional diffusion equatiorPIH 



(4) 



with /3 = 1, so that the time derivative on the left hand side is a first- 
order derivative. The fractional space derivative on the right hand side 
is a Weyl-Rietz fractional derivative^ defined below. Here the anomalous 
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diffusion coefficient has units cm'^/sec. A fundamentai cfialienge is to 
derive sucli fractionai equations from a microscopic tfieory, without invok- 
ing power-law statistics in the first place. Furthermore, the solutions of 
such equations exhibit a diverging mean-square displacement (x^) — oo, 
which violates the principle of causalitjQ which restricts physical phenom- 
ena to spread at finite speeds. So how can fractional equations like Eq. 
Q describe physical reality? We will address this paradox in this work. 
The solution of Eq. Q for an initial narrow cloud is given in terms of a 
Levy distribution (see details below) . The Levy distribution generalizes the 
Gaussian distribution in the mathematical problem of the sum of a large 
number of independent random variables symmetrically distributed about 
0, in the case where the variance of the summands diverges, corresponding 
physically to scale free systems. Here our aim is to derive Levy statistics 
and the fractional diffusion equation from the semi-classical picture of Sisy- 
phus cooling. Specifically, we will show that (3 — 1 and relate the value of 
the exponent v to the depth of the optical lattice Uq, deriving an expres- 
sion for the constant K,^. Furthermore, we discuss the limitations of the 
fractional framework, and show that for a critical value of the depth of the 
optical lattice, the dynamics switches to a non-Levy behavior (i.e. a regime 
where Eq. Q is not valid); instead it is related to Richardson-Obukhov 
diffusion found in turbulence. Thus the semiclassical picture predicts a rich 
phase diagram for the atomistic diffusion process. We will then compare 
the results of this analysis to the experimental findings, and see that there 
are still unresolved discrepancies between the experiment and the theory. 
Reconciling the two thus poses a major challenge for the future. 

Since usual approaches to diffusion, namely Fick's second law and the 
Green-Kubo formula break down for the case of interest, we cannot use 
ordinary approaches. Here the power of first passage time statistics enters. 
We will show how analysis of first passage time statistics yields the diffusion 
front of a packet of particles. In essence the first passage time tool replaces 
the failing Green Kubo formula, and leads us to the solution of the problem, 
namely the exponents /3 and v as well as the generalized diffusion constant 
K^. But this is a long journey, so first let us briefiy review the concept 
of Brownian excursions, first passage times for simple Brownian motion, 
and the distribution of the area swept under Brownian motion till its first 
passage. 



''See the discussion in Ref. |10| where the unphysical nature of Levy flights is discussed 
and the resolution in terms of Levy walks is addressed. 
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1.2. Brief Survey: Joint PDF for first passage time and the 
area swept under Brownian motion 

Consider a Brownian motion in one dimension. The time Tf it takes the 
particle starting at a; = e > to reach a; = for the first time is the first 
passage time. The distribution of t/ has been well investigated both for 
Brownian motion and for other types of random walksISES A quantity 
which until recently was only of mathematical interest is the total area 
under the Brownian path in the time interval (0, r^), which was treated 
by Kearney and Majumdar.^ This area, which we denote Xf: is obviously 
positive if we start at a; = e > 0. Later we will take e — )■ which is a 
subtle point, but for now we keep e finite. The pair of random quantities 
Xf and Tf are clearly correlated since a large Tf implies statistically a large 
Xf- The joint probability density function is denoted as P(Tf,Xf)- Using 
Bayes' law we may write this density as 



where g{Tf) is the first passage time probability density function (PDF), 
which at least for Brownian motion is well investigated. In particular, gijf) 
decays as a power for large Tf, 



so that the average first passage time, (r/) = cx), diverges for one dimen- 
sional Brownian motion on the infinite line. P{Xf\'''f) is the conditional 
probability density function of %/ given some fixed Tf. This leads us to the 
concept of a Brownian excursion. 

A Brownian excursion is a conditioned one dimensional Brownian mo- 
tion x{t) over the time interval < t < (recall here Tf is fixed to a 
specific value). The motion starts at x{{)) = e > 0, e <C 1, and ends at 
x{Tf) — and is constrained not to cross the origin a; = in the observa- 
tion time {0,Tf), namely the Brownian excursion is always positive. The 
area under the Brownian excursion is the object of interest since that gives 
the conditional probability density Pixfhf)- Thus the statistics of the area 
under the Brownian excursion also gives the joint PDF of the first passage 
time and the total area under the Brownian motion. 

The area under the Brownian excursion was treated by mathematicians 
with applications in computer science and graph theory and more recently 
in the physical context of fluctuations of interfaces Majumdar and 
ComtePI give a concise introduction and derivation of the solution to this 
problem in terms of the Airy distribution. We shall soon show that the 



^{Tf,Xf) ^ 9{Tf)P{Xf\Tf) 



(5) 



giTf) cc {Tf)-^'\ 



(6) 
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Fig. 2. An example of a Brownian excursion x{t), which in the text is generaUzed to a 
Langevin excursion in momentum space, p{t'). These are random curves which start and 
end close to the origin and remain positive for a fixed time. Surprisingly the random area 
under these curves has started to attract physical attention in diverse problems ranging 
from fluctuating interfaces and as we show here diffusion of atoms in optical lattices. A 
related problem is the Tracy-Widom distribution. 

solution of the dynamics of the atomic packet demands the solution of a first 
passage time in momentum space, and the distribution of the area under the 
random first passage time curve. Thus the random area under a random 
curve terminated at its first passage turns into a relevant problem, but, 
fortunately or unfortunately, the apposite curve is not a simple Brownian 
one. 

So, for perspective, we first review what is known about the distribution 
of A, the random area under the Brownian excursion, A = J^^ x{t)dt in the 
e limit. Since Brownian motion scales like x ~ f^/^ the area A scales 
like J^^ x{t)dt cx (r/)^/^. Thus the distribut ion has the scaling behavior 
Prj(A) = {Tf)~^/'^ f{A/{TfY/'^), where the (t/)~^/^ prefactor stems from 
the normalization condition. The Laplace x ^ u transform of /(x), f{u) 
was computed by Darling and LouchardHEZl 

oo 
k=l 
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(here, the diffusion constant of the Brownian motion, arising from the sim- 
ple random walk used to model the motion, is 1/2). The a^s are the 
magnitudes of the zeros of the Airy function Ai{z) on the negative real 
axis ai = 2.3381.., a2 = 4.0879..., — 5.5205.. and for that reason the 
distribution is called the Airy distribution. TakacJ^^ inverted the Laplace 
transform to find 

/(^) = ^E^-'^-^^'^f t^("5/6,4/3,6./x^) (8) 
fc=i 

where hk = 2(Q;fe)'^/27 and [/(a, 6, z) is the confluent hypergeometric func- 
tion.liSl Moments and asymptotic behaviors of the Airy distribution can be 
found in Ref. |T5] The function can be calculated and plotted by truncating 
the sum at some x dependent value. 

From this, we can deduce the power law exponent describing the tail of 
the PDF of the area swept under the Brownian curve till its first passage 
event ^(x/), averaged over all Tf. Clearly the marginal PDF q{xf) is found 
from the joint PDF in the usual way by integration 

/>oo 

dixs) = / g{Tf)P{xf\Tf)dTf (9) 

Using Eq. (j6| and the scaling behavior of the conditional PDF P(x/|r/), 
we have 

/>oo 

9(x/)«/ {r})-''^{Ti)-''^f{xfl{rsf'^)dT;. (10) 

Changing integration variables to y = Tf / {xf)'^^^ g^t the asymptotic 
form 

9(X/) « iXf)-'"- (11) 
We see that the scaling exponent 3/2 for g{T) and 4/3 for q{xf) are easy 
to obtain and are both related to the scaling of Brownian motion x cx t^/"^. 
The full distribution of q(x/) was recently found by Kearney and Majumdar 
for finite e. We will later deal with such a calculation in more rigor, but 
for now we just focus on the exponents. We now turn to showing how the 
generalization of these concepts lead to the description of the dynamics of 
cold atoms. 

2. Model and Goal 



We wish to investigate the spatial density of the atoms, P{x, t). The trajec- 
tory of a single particle is x{t) = p{t)dt / m where p{t) is its momentum. 
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Within the standard pictureSH of Sisyphus coohng, two competing mecha- 
nisms describe the dynamics. The coohng force 

F{p) = -ap/[l + {p/p,f] (12) 

acts to restore the momentum to the minimum energy state p = Q. Mo- 
mentum diffusion is governed by a diffusion coefficient which is momentum 
dependent, D{p) = Di + D2/[1 + [p/pcY]. The latter describes momentum 
fluctuations which lead to heating (due to random emission events). We 
use dimensionless units, time t — ta, momentum p — )■ p/pc, the momentum 
diffusion constant D = Di/ (pc)'^a and x — xma/pc- For simphcity, we set 
D2 — since it does not modify the asymptotic |p| — oo behavior of the 
diffusive heating term, nor that of the force, and therefore does not modify 
our main conclusions. The Langevin equations 

^ = F(p) + -mm, (13a) 

describe the dynamics in phase space. Here the noise term is Gaussian, 
has zero mean and is white {S,{t)£,{t')) = 5{t — t'). The now dimensionless 
cooling force is 

^b)-rT^. (14) 



The stochastic Eq. (131 give the trajectories of the standard Kramers 
picture for the semi-classical dynamics in the optical lattice which in turn 
was derived from microscopical considerationsEE! From the semiclassical 
treatment of the interaction of the atoms with the counterpropagating laser 
beams, we have 

D = cEr/Uo, (15) 

where Uo is the depth of the optical potential and Er the recoil energy, and 
the dimensionless parameter c depends on the atomic transition involved]^ 
For p <^1, the cooling force Eq. ( |l4| ) is harmonic, F{p) ~ — p, while in the 
opposite limit, p ^ 1, F{p) ^ — 1/p. It is useful, given the correspondence 



between Eq. (13a) and overdamped motion in space subject to a force F, 



''For atoms in molasses with a. Jg = 1/2 — » Je = 3/2 Zeeman substructure in a lin 



. 1 /2 

lin laser configuration c = 12.3.1^ Refs. |5|20|2l| give c = 22. As pointed out irP' different 
notations are used in the literature. 
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to define an effective potential 

V{p) = - r F{p)dp = (1/2) ln(l + p2). (16) 
Jo 

This potential is asymptotically logarithmic, V{p) ~ ln(p) when p is 
large, which has as its consequence various unusual equilibrium and non- 
equilibrium properties of the momentum distributionP^HSU pgr example the 
steady state momentum distribution is Weq{p) oc exp[—V{p)/D] which has 
a power law form 

W,gip)=Af{l+p')-'^^° (17) 

where Af is the normalization. Such non-Gaussian behavior was observed 
in the experiments of Renzoni's group. For D < 1/3, corresponding to 
shallow optical lattices, the averaged kinetic energy diverges. A similar 
effect was found experimentally by Walther's groupP^ when the depth 
of the optical lattice is tuned one may observe a gigantic increase of the 
energy of the system. Of course the energy of a physical system cannot be 
infinite and this paradox was recently resolved ? '^ I ^^ l by treating the problem 
dynamically, revealing that the energy may increase with time but it is never 
infinite. Of course, when the energy becomes too large, other experimental 
effects become relevant as well. 

While the velocity distribution of the atomic cloud appears to be well un- 
derstood, the new experimenlP demands a theory for the spatial spreading. 
This will be the focus of this chapter. However, by now the reader might 
think that laser cooling is an oxymoron. Namely if we get non-Maxwellian 
distributions in equilibrium, fat tails, and Levy behavior, where is the ther- 
mal state found? The answer, as we understand it, is rather involved. First 
the equilibrium momentum PDF can be rather narrow in its center, and in 
that sense the atoms are cold. Further, the divergence of energy happens 
for low D, so an experimentalist interested in cooling will turn their knobs 
away from that limit. In fact there is an optimal depth of the optical po- 
tential in the sense that it yields the smallest value for the mean kinetic 
energy of the particles, since the basic energy scale is set by the optical 



potential depth (see, e.g.. Fig. 4 of Ref. j22jand Fig. 2 of Ref. |26|). Strictly 
speaking, however, for all D the equilibrium state is not thermal since the 
momentum distribution is not Gaussian. Experimentalist intuitively rec- 
ognize this problem, and hence also use a method called evaporation to 
reach an ultimate cold state. Evaporation simply means that atoms in the 
tails, i.e. atoms with high velocities, are a problem for cooling, so they are 
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allowed by the experimentalist to leave the system. Another mechanism 
that may lead to a realistic thermal state is internal collisions among the 
atoms (which will depend of-course on the density). It is believed that 
such effects will lead the system in the long time limit to thermal equi- 
librium. However, at least in the experiments briefly surveyed herePEESl 
which take long times on atomic scale (e.g. milliseconds) no hint of thermal 
Maxwellian equilibrium is found, at least for values of D of interest. Hence 
it is an important task to map the different phases of the dynamics which 
are controlled by D. Indeed this system is truly wonderful in the sense that 
the experimentalist can tune D and hence explore different effects both in 
and out of equilibrium. 




Fig. 3. Schematic presentation of the momentum of the particle versus time. The 
times between consecutive zero crossings are called the jump durations r and the shaded 
area under each excursion are the random flight displacements x- The t's and the x's 
are correlated, since statistically a long jump duration implies a large displacement. 



3. From Area under the Bessel excursion to Levy walks for 
the packet of atoms 

The heart of our analysis is the mapping of the Langevin dynamics to a 
recurrent set of random walks. The particle along its stochastic path in 
momentum space crosses p — many times when the measurement time 
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is long. Let r > be the random time between one crossing event to the 
next crossing event, and let — oo < x < oo be the random displacement 
(for the corresponding r). As schematically shown in Fig. |3j the process 
starting at the origin with zero momentum is defined by the sequence of 
jump durations, {ti, T2, . . .} with corresponding displacements {xi, X2, • ■ •}, 
with xi = Jq^ p{T)dT, X2 = Jri^'^^ P(^)'^'^^ total displacement x 

at time i is a sum of the individual displacements Xi- Since the underlying 
Langevin dynamics is continuous, we need a more precise definition of this 
process. Let t be the time it takes the particle with initial momentum pi to 
reach pf ~ Q for the first time. So r is the first passage time for the process 
in momentum space. Eventually we will take pi ~> Pf, see below. Similarly, 
X is the displacement of the particle during this flight. Namely by definition 
X is the area under the Langevin curve (in momentum space) till the first 
passage time. The probability density function (PDF) of the displacement 
X is denoted g(x) and of the jump durations g(r). To conclude we see 
that first passage time statistics, in particular the area swept under the 
random momentum process p(r) until the first passage time, constitute the 
microscopical random jumps and correlated random waiting times which 
eventually give the full displacement of the the particle. 

As shown by Marksteiner, et al. and LutzP22l these PDFs exhibit power 
law behavior 

g (r) cx r~5~2fe, (18a) 

9(x)oc|xr^"^, (18b) 

as a consequence of the logarithmic potential, which makes the diffusion for 
large enough p only weakly bounded. It is this power-law behavior, with 
its divergent second moment of the displacement x for D > 1/5, which 
gives rise to the anomalous statistics for x (measured in the experiment). 
Notice that taking the limit Z? — 00 in Eq. ( [T8| we get the expected 
results as given in Eqs. ^ and (11). In this limit, the cooling force F{p) 



is negligible, the dynamics is governed by the noise in Eq. ( 13 ) and so the 
process p{t) approaches a Brownian motion. This corresponds to atoms 
emitting a photon in random directions, and hence from the recoil jolts, 
they are undergoing an unbiased random walk in momentum space (which 
is the pure heating limit). 

Importantly, and previously overlooked, there is a strong correlation be- 
tween the jump duration t and the spatial extent of the jumps x- These 
correlations are responsible, in particular, for the finiteness of the moments 
of P{x,t). Physically, such a correlation is obvious, since long jump dura- 
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tions involve large momenta, which in turn induce a large spatial displace- 
ment. The theoretical development starts then from the quantity ■'/'(Xi''")) 
the joint probability density of \x\ and r. From this, we construct a Levy 
walk schem^H22l which relates the microscopic information V'(X;''') to the 
atomic packet P{x, t) for large x and t. 

4. Scaling Theory for Anomalous DifTusion 

We rewrite the joint PDF 

V'(x,r) = ^5(r)p(|xl|T), (19) 

where p{x \ t) is the conditional probability to find a jump length of < 
X < oo for a given jump duration r. The factor ^ accounts for the walks 
with X < 0. Numerically we observed that the conditional probability scales 
at large times like 

p(x|T)«2T-^S(|x|/r^); 7 = 3/2. (20) 

B{-) is a scaling function soon to be determined. To analytically obtain the 
scaling exponent 7 = 3/2 note that q (x) = /q°° dr?/; (Xi ''")j giving 

'?(x)^^"dTr-i-riTr-^i?(^M^ cx|xr(^+'^). (21) 
Here tq is a time scale after which the long time limit in Eqs. (|18[) holds 



and is irrelevant for large x- Comparing Eq. (21) to Eq. ( |18b[ ) yields 
the consistency condition 1 + (1 + l/D)/{2'^) = 4/3 + 1/{3D) and hence 
7 = 3/2, as we observe in Fig. |4] 



We see that the scaling function B{.) in Eq. (20) is similar to the 
problem of the area under a Brownian excursion. However, now we are 
confronted with the additional impact of the friction force. The concept of 
the Brownian excursion needs to be extended to a stochastic curve called 
the Bessel excursioUjI^Sl which corresponds to with the friction force F(j)) = 
— 1/p. The term Bessel excursion stems from the fact that mathematically 
speaking, diffusion in momentum space, in the non-regularized potential 
ln(p), corresponds to a process called the Bessel process l^il22J As with the 
Brownian excursion, the Bessel excursion is the Langevin path p{f) over 
the time interval < t' < t such that the path starts on — > and ends 
on the origin, but is constrained to stay positive (if pi > 0) or negative 
(if Pi < 0). Brownian excursions are just the D oo limit of this more 
general problem. Schematic excursions p{t) are presented in Fig. [sj albeit 
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Fig. 4. Lower panel: We plot the flight displacement |x| versus the jump duration r to 
demonstrate the strong correlations between these two random variables. Here D = 2/3 
and the line has slope 3/2 reflecting the x ^ t^/^ scaling discussed in the t ext. U pper 
panel: The conditional probability p{x\t) is described by the Airy distributiorEElwhen 
£> — >■ oo; otherwise by the distribution derived here based on the Feynman-Kac theory, 
Eq. l |24| . They describe areas under constrained Brownian (or Bessel) excursions when 
the cooling force F(p) = (or F{p) ^ 0) respectively. 



there with various excursion times r. The displacement x = Jq P{^'')'^t' is 
the area under the excursion as shown in Fig. [Sj 



May 2, 2013 2:10 World Scientific Review Volume - 9in x 6in ws-Barl<aiKes5ler2 



Transport and the first passage time problem with application to cold atoms in optical trapsl5 

4.1. Area Under the Bessel Excursion 

We briefly outline how to obtain the distribution of the area under the 



Bessel excursion, namely we find the rather formidable p{x\''')^ Eq. (24) 
below. The main tool we use is a modified Feynman-Kac theory. The origi- 
nal well-known formalism deals with functionals of Brownian motion, while 
here we deal with functionals of Langevin motion, hence the modification. 
The Feynman-Kac equation is essentially a Shrodinger equation in real time 
and is reviewed in Ref. |33[ The calculation involves three steps: 

(a) We first find G{s,p, r), the Laplace transform with respect to x of 
the joint PDF for the positively constrained path to reach p with 
area x time r, which we denote by G{x,P, t). 

(b) Taking the limit of the initial and final momentum pi ^ p ^ 0, we 
get p{s\t) = limp^p.^o G'(s,p,r)/G(s = 0,p, r); the G{s = 0,p, r) 
denominator yields the correct normalization. 

(c) Finally, using the inverse Laplace transform, we get p{x 1 1). 

As mentioned, the main tool we use is the modified Feynman-Kac formalism 
for functionals of over-damped Langevin pathsEHMl For step (a) we solve: 



dGjs.p, t) 
dr 



Lfp-spJ G(s,p,t), (22) 
where Lfp = D{dpY — dpF{p) is the Fokker Planck operator. When F{p) = 



0, Eq. (22 1 is the celebrated Feynman-Kac equation which is an imaginary- 
time Schrodinger equation in a linear potential (i.e. the —sp term). We do 
not provide here a derivation of Eq. (33) but in Sec. [s] provide a derivation 



of a related equation. The constraint p > gives the boundary condition 
G{s,p — 0,t) = 0, i.e., an infinite potential barrier when p < in the 
quantum language. We use the force F{p) = —1/p since we are interested 
only in the scaling behavior of the problem where large excursions imply 
that the small p behavior of the force is not important. This assertion can 
be mathematically justified, but here for the sake of space we will only 
demonstrate it numerically (see Fig. |4]) . Solving Eq. ( 22 ) and following 
the recipe (a) - (b) we find 



fc=i ^ ^ 



(23) 
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Here Dg'l + (gk/pY ~P9k = —DXkgu is the eigenvalue equation which gives 
the eigenvalues and eigenstates gk{p) which are normalized according 
/o°°[fffe(p)]^P^^^'^P = 1. A detailed mathematical derivation will soon be 
published!^ 

It is easy to tabulate the eigenvalues and the slopes on the origin 
d'kiv ~ 0) with standard numerically exact techniques. With this infor- 
mation and the inverse Laplace transform we find the solution in terms of 
generalized hypergeometric functions 



p{x 



r(i 
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(24) 



where 



273 



2.7 



I' sm TT 



3 



J=0. 



(25) 

Here we have introduced the 



and the summation is over the eigenvalues 
parameter 

l + D 

V — 

31) 

which will turn out to play a crucial role in the dynamics 
X ~ T'^/2 scaling which proves the scaling hypothesis, Eq. ( 20 1 



(26) 

Notice the 
To compare 

with numerical simulations, we evaluate this function truncating the sum 
beyond some k, since for fixed x the terms decay rapidly in k. The Langevin 
simulations were performed by mapping the problem onto a biased random 
walk in momentum space, with momentum spacing 5p = 0.1. 

In Fig. |4j we present p{x\t) obtained from numerical simulations, show- 
ing that it perfectly matches our theory on Bessel excursions without fitting. 
Our theory and simulations show that the areas under the Brownian and 
Bessel excursions share the same x ~ t^^"^ scaling behavior. In both cases, 
p(x|t) falls off rapidly for large |x|, ensuring finite moments of this distri- 
bution, which will prove important later. When ^ 1, the Bessel and the 
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Brownian excursions coincide, since then the force F(p) plays a vanishing 



role, as can be shown directly from Eq. (24). The Bessel excursion PDF for 
finite D is compared with the limiting Airy distribution in Fig. |4j where it 
is seen that the Bessel excursions are pushed further away from the origin. 
The —1/p force implies that trajectories not crossing the origin in (0, t) are 
typically further away from the origin, compared with the free Brownian 
particle, since trajectories close to the origin are more likely to cross it, due 
to the attractive force. 

4.2. Coupled Continuous Time Random Walk Theory from 
First Passage Time Statistics 

The next step is to relate the first passage time properties of the under- 
lying Langevin process in momentum space, namely the joint distribution 
of the area under the Bessel excursion and the first passage time, with the 
spreading of the density P{x,t). Given our scaling solution for p{x\t), and 



^(t) which is considered in detail in Ref. 35 (see also Appendix here) 
and hence ip{x,t), we construct a theory for P{x,t) using tools devel- 
oped in the random walk community.'^^l As soon detailed, one first obtains 
a Montroll-Weis^ type of equation for the Fourier-Laplace transform of 
P{x,t), P{k,u), in terms of tp{k,u), the Fourier-Laplace transform of the 
joint PDF V(x,t): 

Here, '9{k,u) is the Fourier-Laplace transform of t^^^'^B{\x\/t'^^^)[1 — 



/p ■!/'(r)(ir]. The last step is then to invert Eq. (271 back to the x, t 
domain. 



To derive Eq. (27) we present a formalism which relates the joint dis- 
tribution 

V'(x,r)=5(r)r-3/2B(|x|/T3/2) (28) 
with P{x,t) {B{.) is given in Eqs. ( 20|24 )). Define r]s{x,t)dtdx as the 



probability that the particle crossed the momentum state p = for the 
s"' time in the time interval (t, t + dt) and that the particle's position was 
(x, x+dx) . This probability is related to the probability of the s— 1 crossing 
according to 



Vsix,t) = J dv3/2 J dTr]s-i(x ~V3/2T^^'^,t- B (\v3/2\) dir) 

(29) 
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where we changed the variable of integration to = xl^'^^'^ ■ Now the 
process is described by a sequence of jump durations ti, T2, • • • and the cor- 
responding generalized velocities 'y3/2(l), ^3/2(2), • • • (with units [a;]/[t^/^]). 
The displacements in the sth interval is: Xs — ''^3/2(s)[ts]'^^^- The advan- 
tage of the representation of the problem in terms of the pair of microscopic 
stochastic variables r, W3/2 (instead of the correlated pair r, x) is clear from 



Eq. (29 1: we may treat 113/2 and r as independent random variables whose 
corresponding PDFs are gir) and Biy-^j-i) respectively. The initial condi- 
tion a; = at time t = implies r]o{x,t) = S{x)S{t). Let P{x,t) be the 
probability of finding the particle in {x, x + Ax) at time i, which is found 
according to 

/"OO 1*00 

P{x,t)^y2 ^^3/2/ dTfjJx~Vy2T^/^,t-T)B{\vy2\)W{T). 

(30) 

since for a time series with s intervals, the last jump event took place 
at i — r and in the time period (t — r, t) the particle did not cross the 
momentum origin (hence the survival probability W{t) — 1 — Jj^5(r)dr). 



The summation in Eq. ( 30 1 is a sum of all possible realizations with s 
returns to the momentum origin pf = 0. 

As usual we consider the problem in Laplace-Fourier space where 



t ^ u and X ^ k. Eq. (29) then translates to 

77s(fc, u) = 77s-i(fc, u) '4j{k, u). (31) 
Here the tildes refer to the Fourier-Laplace transforms, and we have used 



Eq. (20) to write that ipiXiT) — gi''')''' ^^'^B{\x\/t'^^^)- Hence, summing 



the Fourier-Laplace transform of Eq. (30) using Eq. (|3l|), we get the 



MontroU- Weiss Eq. (27 1 



The Montroll- Weiss equation is widely used in applications of the ran- 
dom walk. This equation holds since the underlying dynamics is Markovian, 
namely each time the process p{t) crosses zero, the process is renewed. How- 
ever, there is a subtle point behind our approach, which makes it different. 
Since the trajectories are continuous, the number of zero crossings s is in- 
finite within any finite time interval (0,t), and for a particle starting with 
p — 0. In contrast in usual continuous time random walk theories the num- 
ber of renewals is always finite, and it increases with time. The functions 



entering the right hand side of Eq. (27 1 still depend on the initial momen- 



tum we assign to the particle ±pi immediately after each zero crossing (if 
we use a grid in a numerical simulation this is handled automatically) . We 
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still must take the limit Pi — > and show that our final expression will give 



sensible results in that regularized limit (see also Ref. 35 1 . 



5. The Levy phase 

We now explain why Levy statistics describes the diffusion profile P{x,t) 
when 1/5 < D < 1, provided that x is not too large. The key idea is that, 
for x's which are large, but not extremely large, the problem decouples, and 
^/'(fc, u) can be expressed as a product of the Fourier transform of q(x), q{k) 
and the Laplace transform of ^(t), g{u). This is valid as long as x ^ t^^^, 
since otherwise paths where x ~ t^^^ are relevant, for which the correlations 
are strong, as we have seen. The long-time, large- a; behavior of P{x, t) in 
the decoupled regime is then governed by the small-fc behavior of q{k) and 
the small u behavior of g{u). When the second moment of q{x) diverges, 
i.e. for D > 1/5, the small-A; behavior of q{k) is determined by the large-x 



asymptotics of q{x) as given in Eq. (20), q{x) ~ x'^ /\x\^~^'^ ■ When the 
first moment of r is finite, i.e. for D < 1, the small-it behavior of g{u) is 
governed by the first moment, (r). From these follow the small-fc, small-u 
behavior of P(fc, u): 

where K^, = 7ra;^/((T)r(l + ly) sin ^) (see some details in Appendix A). 

Both x'^ and (r) can be calculated (see Appendix A and Ref. |35j) via 
appropriate backward Fokker-Planck equations. They both vanish as the 
magnitude of the initial momentum of the walk goes to zero, but their ratio 
has a finite limit, so that K^, upon returning to dimensionfull units, is 



r(^)32— i[r(z/)]2sin(^) 



m 



P{k, u), as given in Eq. ( |32| , is in fact precisely the symmetric Levy distri- 
bution in Laplace-Fourier space with index i', whose (x, t) representation is 



(see Eq. (B17) of Ref. 36]) 



P{x,t) 



(34) 



The properties of the Levy function L^fl{.) are well-known, and the solution 
in time for the spatial Fourier transform is P{k,t) — ex-p{—K^t\k\'^) with 
2/3 < < 2. This distribution is the solution of the fractional diffusion 
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equation, Eq. Q, with (i — 1 and an initial distribution located at the 
origin. This justifies the use of Eq. Q in Ref. 4 for 1/5 < Z? < 1 and 
provides v and Ki, in terms of the experimental parameters. We can verify 
this behavior in simulations, as shown in Fig. 3, where we see excellent 



agreement to our theoretical prediction, Eqs. (34) and (33), without any 
fitting. 

The fractional diffusion equation, Eq. [4] has the following meaning 



and connection to the underlying random walk. We re-express Eq. ( 32 ) 
as uP{k,u) - 1 = -K^lkl^Pik.u). Recalling that the Fourier space rep- 
resentation of the Weyl-Rietz fractional derivative, V^, is — |fc|'' and that 
uP{k,u) — 1 is the representation of the time derivative in Laplace space 
for i5-function initial conditions centered at the origin, we see that in fact 
P{x,t) satisfies the fractional diffusion equation, Eq. (|4|, with (3 — 1. 

Fig. [5] (b) also illustrates the cutoff on the Levy distribution, which is 
found at distances x ~ t^/'^ , as can be seen also from Fig[5|c). Beyond this 
length scale, the density falls off rapidly (Fig. [5jc)). This, as noted above, 
is the result of the correlation between x ^nd r, as there are essentially no 
walks with a displacement greater than of order t^^"^. This cutoff ensures 
the finiteness of the mean square displacement, using the power law tail of 
the Levy PDF L^{x) ^ x-'^^+''^ and the cutoff we get: 



/ t-(i/-)(a;Ai/'')-(i+'^)x2da;^t4-3./2^ 



for 2/3 < < 2, in agreement with rigorous results of Ref. 37 (where the 
prefactor is also computed). As noted in the introduction, if we rely on the 
fractional diffusion equation, Eq. naively, we get (x^) — 00. Thus the 
fractional equation must be used with care, realizing its limitations in the 
statistical description of the moments of the distribution and its tails. 

6. Richardson Phase 

When the average jump duration, (r), diverges, i.e., for D > 1, the dynam- 
ics of P(x,t) enters a new phase. Since the Levy index v approaches 2/3 
as D approaches 1, x scales like t^/"^ in the limit. Due to the correlations, 
X cannot grow faster than this, so in this regime, P{x,t) ^ t~^/'^h{x/t^/'^). 
This scaling is that of free diffusion, namely momentum diffuses likep ~ t^^'^ 
and hence the time integral over the momentum scales like x ^ t^^'^. In- 
deed, in the absence of the logarithmic potential, namely in the limit Z? 3> 1 
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Fig. 5. (a) tl/'"P(|x|,t) versus \x\/t^^'' for D = 2/5, i.e. = 7/6. The theory, i.e., 
the Levy PDF Eq. ( |34[ | along with K^, Eq. ( |33[ l, perfectly matches simulations without 
fitting, (b) Same data in log-log scale. Notice the cutoff for large x which is due to 
the coupling between jump lengths and jump durations, making x > t^^^ extremely 
improbable, (c) t^^'^^^')^^P(\x\,t) versus |a;|/t''/^ for D = 2/5, showing the crossover 
from power-law to Gaussian behavior at \x\ ~ t^^^ . This cutoff ensures the finiteness of 
the mean square displacement of the underlying spreading of the packet of atoms. 



Eq. (13) gives 



exp 



3a; 



2 1 



(35) 



This limit describes the Obukhov model for a tracer particle path in tur- 
bulence, where the velocity follows a simple Brownian motionlM2£l These 
scaling properties are related to Kolmogorov's theory of 1941 (see Eq. (3) 
in Ref. |39| and to Richardson diffusion^^ ^o] ^^2^ ^ ^3 ^ jg yg^jjj 

when the optical potential depth is small since Z? — > oo when Uq — > 0. This 
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limit should be taken with care, as the observation time must be made large 
before considering the limit of weak potential. In the opposite scenario, i.e. 
Uq ^ Q before t — >■ oo, we expect ballistic motion, |x| ^ t, since then the 
optical lattice has not had time to make itself felt.^^ 

7. Relation between first passage times and diffusivity 

When the variance of x is finite, namely u > 2, we get normal diffusion so 
that P{x,t) is a Gaussian. In that case the spatio-temporal distribution 
of jump times and jump lengths decouples ipiXiT) — q[x)9{T)- The mean 
square displacement (x^) ^ 2K2t and the diffusion constant is 

This equation relates the first passage time properties of the momentum 
process and the spatial diffusivity and it has the structure of the famous 
Einstein relation, relating the variance of jump sizes and the average time 
to make a jump, with a diffusion constant. However here we are dealing 
with a special limit, namely (r) is the average time it takes for a particle 
with initial momentum pi = ±e to reach p = for the first time. This 
average first passage time vanishes in the limit of e — > 0. The variance of x, 
is (x^) — J^oo X^9(x)dx namely the variance of the area under the random 
first passage process p{t') till its first return to p = 0, and it too vanishes 
in the limit e 0. Notice that here x may be either positive or negative 
(see Fig. |3]) and its mean is zero (since the particle may start on pi — ±e) 
with equal probability since the force F{p) is zero on p = 0. Of-course from 
symmetry F{p) — —F(—p), so it suffices to consider x for say positive areas 
only. 

The calculation of (r) and (x^) is based on backward Fokker-Planck 
equations Eqs. ( A.1|A.3 ) (see further details in Ref. 35) and we find 



-1 /"OO r /"OO 



(37) 



and Z = J^^exp[—V{p)/D]dp. This equation was derived previously using 



a different approach in Ref. 41 The coefficient if 2 diverges as D — >■ 1/5 
from below indicating the transition to the superdiffusive phase. A sharp 
increase in the diffusion coefficient of an atomic cloud as the laser intensity 
reaches a critical value was detected experimentally by Hodapp et al.,1^ 
which could be the fingerprint of the transition from Gaussian to the Levy 
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phase. We see that the Green-Kubo formula, Eq. ([3|, is not the only 
way to calculate the diffusion constant. For Markovian processes, with 



non-linear friction, using Eq. (36) and calculating the first passage time 
properties is an alternative. In the normal diffusion case, the Green-Kubo 
formalism is so successful that the existence of an alternate method is only 
of academic interest. However, as shown throughout this current work, 
calculations of first passage time properties of the velocity process are very 
useful for the unraveling of the anomalous spatial diffusion properties of the 
system, where we cannot use the standard Green-Kubo formalism nor the 
Gaussian central limit theorem. See further work on a generalized Green- 
Kubo formalism and non-stationary momentum correlation functions in 
Ref. M 

8. Equation for the distribution of the area swept by a par- 
ticle until its first passage 

We briefly outline the derivation of the equation for the PDF of the area 
X = /p p(i')dt' under the stochastic process p{t) > which starts on pi > 
until its reaches p = for the first time, so that t is the first passage time. 
The PDF (?(|xl) has a Laplace x — >■ s transform which we call Q{s,pi) (here 
X > the extension to negative x being trivial). The equation for this 
quantity is 



D^+F{p^^-sp,Q = 0. (38) 
dpt dpi 



For simplicity and for conciseness we consider the case F{p) = 0, so p{t) is 



undergoing Brownian motion without bias. By definition Xpi = Jq p{t')dt' 



so 



Q{s,p,) = (exp[-s / p(t')dt'])p. (39) 







where t is the first passage time. The starting point is denoted by the 
subscript. Now we consider the process p{t) as a random walk on a lattice 
with probability 1/2 to jump from p(t) — >■ p{t) + 6p or p(t) — > p{t) — Sp, 
so 5p is the lattice spacing. The constant times between these jumps is St 
and we consider the limit Sp ^ and St ^ while keeping D — [Sp)"^ /2St 
finite, this being the Brownian limit. 

A particle starting with pi at time t — Q can either have momentum 
Pi — Sp at time St with probability one half, or pi -\- Sp with the same 
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probability. From the definition of Xpi for a particle starting with pi 
we have 

^ { PrSt + XjH-Sp Prob 1/2 , , 

\p^5t + xp^+sp Prob 1/2. ^ > 

The first case corresponds to a particle initiaUy jumping down in p, pi — >■ 
Pi — Sp while the second case is for a particle jumping up. The term piSt 
stems from the displacement during the interval St, and the remaining term 
(either Xpi or Xpi+25p is the displacement for a new process which starts 
either with pi — Sp or pi + Sp) . Eq. ( |40[ ) is now used to find 

(exp(-sxpj) = ;5 X! (^^Ph^iP^^* - ^Xp^+>jSp])- (41) 
^ ^=±1 



Expanding the right hand side of Eq. (41) to second order in Sp and first 
order in St, and inserting these expansions in Eq. (41 ), using D = (Sp)^ /2St, 



we get Eq. (38). For F{pi) ^ a similar approach is used though now 
weights on the jumping probabilities must be implemented (see similar 
expansions in Refs. 13|34 ). The boundary condition Q{s,Pi)\p-=o = 1 is 



used since a particle starting on — will immediately reach p = and 
then the PDF of x is a delta function centered on zero, so the Laplace 



transform of this delta function Q(s,0) is equal to unity. We use eq. (38) 
in Appendix A (and in more detail in Ref. ^35) to derive the long tail of 
q{x) ^rid from it K^. 

9. Discussion of the experiment 

Our work shows a rich phase diagram of the dynamics, with two transi- 
tion points. For deep wells, D < 1/5, the diffusion is Gaussian, while 
for 1/5 < Z) < 1 we have Levy statistics, and for D > 1 the Richardson- 
Obukhov scaling x ^ t^/"^ is found. So far, experiments detected a Gaussian 
phase and an anomalous one. A transition to Richardson-Obukhov scaling 
has not yet been detected. Sagi et al. ^ clearly demonstrated that changing 
the optical potential depth modifies the anomalous exponents in the Levy 
spreading packet. However, the experiment showed at most ballistic be- 
havior, with the spreading exponent S, defined hy x ^ t^ , always less than 
unity. This might be related to our observation that to go beyond ballistic 
motion, 5 > 1, one must take the measurement time to be very long. A 
more serious problem is that, in the experimental fitting of the diffusion 
front to the Levy propagator, an additional exponent was introducecff to 
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describe the time dependence of the full width at half maximum. In con- 
trast, our semi-classical theory shows that a single exponent v is needed 
within the Levy scaling regime 1/5 < D < 1, with the spreading exponent 
6 = 1/v. This might be related to the cutoff of the tails of Levy PDF 
which demands that the fitting be performed in the central part of the 
atomic cloud. On the other hand we cannot rule out other physical effects 
not included in the semiclassical model. For example it would be very in- 
teresting to simulate the system with quantum Monte Carlo simulations, 
though we note that these are not trivial in the ^ t^/^ regime since the 
usual simulation procedure introduces a cutoff on momentum, which may 
give rise to an artificial ballistic motion. Further, the experiment suffers 
from leakage of particles, which could be crucial, while our theory assumes 
conservation of the number of particles. Another effect our theory does not 
address is the Doppler cooling which may serve as a cutoff at long time 
scales (see some discussion on this in^^). Thus while there is some tanta- 
lizing points of contact between the theory and experiment, full agreement 
is yet far away. 

10. Summary 

Usually first passage times are derived from a random diffusive processes. 
Here our goal was a bit different. We use first passage time statistics to 
derive the spatial diffusivity of the atomic cloud. We revealed also general 
relations between first passage times and the normal diffusion constant K2. 
Since the latter, through linear response, is in many cases related to the 
mobility, we see that analysis of first passage time problems in momentum 
space can be very helpful. For this analysis, we need not only distributions 
of first passage times (which was a topic thoroughly investigated previously) 
but also relatively new tools (at least in physics), namely the distributions 
of areas under the Langevin curve till the random first arrival, and the 
analysis of the area under positive excursion curves like the Bessel excursion 
with fixed time interval. The area under the Brownian excursion, e.g. 
the Airy distribution, is found to be useful in a special limit {D — )> oo) 
which is encouraging since this mathematical object had very little to do 
with physical behavior up until the recent works in the field. The rich 
phase diagram, i.e. Normal, Levy, Richardson diffusions, together with the 
divergence of the diffusion constant close to the transition, revealed by the 
semi-classical approach used here, point to interesting physics, though many 
open questions and challenges both experimental and theoretical remain 
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open with respect to the dynamics and thermodynamics of laser cooled 
atoms. 
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Appendix A. 

In this appendix we derive the small fc, small u behavior of ij:. We first 
not^^H that for v?/"^ ^fc^l,_D<l, we can approximate ip by its 
decoupled form, ip{k,u) ^ g{u)q{k). For small u, since i? is a probability 
distribution, we have g{u) ^ 1 — u(t), as long as the first moment exists. 
This first moment satisfies the backward Fokker-Planck equatiorP^ 

•'^ + ^<'-)? 

dpi dpi 
and the small pi limit of the solution is: 



e 



, - dp^ — A \J (A.2) 

For small k behavior of q{k), given that the second moment diverges, we 
have q{k) ^ 1 — 7ra;^|fc|'^/(r(l + z/)sin^). The length scale x* governing 
the decay of q{x) can in principle be determined from the scaling function 
5, but it is easier to get it directly from another backward Fokker-Planck 
equation for the Laplace transform of q{x)i Q{s,Pi)i where we have made 
explicit the dependence on the initial momentum pi > (so % > 0): 

dpt dp. 

From this, we derive the small pi, small s behavior: Q{s,pi) ~ 1 — piT(l — 
i'){3iy - iys''/{r{i^)3'^''-^), which in turn yields 

<-^- 2.W-ir(.) (^-^^ 

where the 1/2 prefactor stems from equal consideration of positive and 
negative excursions. Putting these two results together, we have 

V^(fc,w)«l-(r)w-<|fc|''--- . (A.5) 



D-^+F{p,)-j^-sp,Q^O . (A.3) 
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To leading order, we may approximate 'i'{k,u) ^ 5'(0,u) ^ ^/{t), which, 



together with Eq. (A.5), gives us Eq. p2 
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